State mapping and discontinuous entanglement transfer in a multipartite open system 
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We describe the transfer of quantum information and correlations from an entangled tripartite 
bosonic system to three separate qubits through their local environments also in the presence of 
various dissipative effects. Optimal state mapping and entanglement transfer are shown in the 
framework of optical cavity quantum electrodynamics involving qubit-like radiation states and two- 
level atoms via the mediation of cavity modes. For an input GHZ state mixed with white noise we 
show the occurrence of sudden death and birth of entanglement, that is discontinuously exchanged 
among the tripartite subsystems. 

PACS numbers: 03.67.Mn,42.50.Pq 



I. INTRODUCTION 

As early as in 1935 Einstein, Podolski, and Rosen [l| 
as well as Schrodinger [2| drew the attention on the cor- 
relations in quantum composite systems and the prob- 
lems raised by their properties. Much later, theoretical 
3 and experimental [|| cornerstones elucidated the is- 
sue of nonlocality. Entanglement is currently viewed as 
the key resource for quantum information (QI) process- 
ing ||, where it allowed a number of achievements such 
as teleportation [6j , cryptography [7j and enhanced mea- 
surements 0. The deep meaning of multipartite entan- 
glement, its quantification and detection [9j, the possible 
applications, are the object of massive investigation. As 
a matter of fact, optical systems have been a privileged 
framework for encoding and manipulating quantum in- 
formation, since bipartite and multipartite entanglement 
may be effectively generated either in the discrete or con- 
tinuous variable regime. On the other hand, the develop- 
ment of QI also requires localized registers, e.g. for the 
storage of entanglement in quantum memories. 

Cavity quantum electrodynamics (CQED) [lfj is a rel- 
evant scenario for this kind of investigations and has been 
addressed for the analysis of entanglement transfer of bi- 
partite El |ll pi 02 [H, [H, 03, El and multipar- 
tite 2(J, 2lL I22LI23]] entanglement. In this framework we 
present a complete study on the entanglement dynam- 
ics of a nine parties open system whose implementation 
could be feasible in the optical regime of CQED [H, [25[ . 
In particular we describe a system where three radiation 
modes, prepared in qubit-like entangled states, are cou- 
pled by optical fibers to three separate optical cavities 
each of them containing a trapped two-level atom. This 
paradigmatic example allows us investigating multipar- 
tite entanglement transfer and swapping in a more real- 
istic way than in [2 CI [2l| , shedding light on fundamental 
processes related to quantum interfaces and memories in 
quantum networks (23 . [231 ] . 

We demonstrate that a complete mapping of pure en- 
tangled states onto the tripartite atomic subsystem oc- 
curs when the external field has fed the cavities. If this 



field is then switched off, the quantum correlations can 
be periodically mapped onto the tripartite atomic and 
cavity mode subsystems according to a triple Jaynes- 
Cummings (JC) dynamics 26]. In the case of exter- 
nal radiation prepared in a mixed Werner state we deal 
with the recently observed phenomenon of entanglement 
sudden death (ESD) (and birth (ESB)) [III, S3, that in 
our case involves the abrupt vanishing (and raising) of 
quantum correlations in tripartite systems. Though this 
is in general a still open problem, nevertheless we can 
show the occurrence of discontinuous exchange of tripar- 
tite entanglement via ESD effects. We also describe the 
dissipative effects introduced by the presence of external 
environments, such as the decay of cavity modes, atomic 
excitations, and fiber modes. Some results are then re- 
ported in the case that the coupling between external 
and cavity mode fields is generalized from monomode to 
multimode. 

In Sec. II we introduce the model of the physical sys- 
tem. In Sec. Ill we derive all main results concerning 
state mapping, entanglement transfer, and entanglement 
sudden death, also adding the presence of external en- 
vironments. The case of multimode external-cavity field 
coupling is addressed in Sec. IV. Some conclusive re- 
marks are reported in Sec. V. 



II. MODEL OF THE PHYSICAL SYSTEM 

We consider an entangled tripartite bosonic system (f), 
prepared in general in a mixed state, interacting with 
three qubits (a) through their local environments (c) also 
in the presence of dissipative effects. In the interaction 
picture the system Hamiltonian has the form: 

J=A,B,C 

+ E hMWjfh + ^jfK)}} (i) 

J,K=A,B,C 



2 



The operators cj, Cj (/j, /j) are the annihilation and cre- 
ation operators for the cavity (external radiation) modes, 
while o-j,<7j are the raising and lowering operators for 
the atomic qubits in each subsystem (J =A,B,C). We 
consider real coupling constant gj for the atom-cavity 
mode interaction and vj,i((t) for the interaction of each 
cavity mode with the three modes of driving external ra- 
diation. We take time dependent constants in order to 
simulate the interaction switching-off at a suitable time 
T ff for the external field. Now we take into account 
three processes of dissipation: the cavity losses at rate 
k c due to interaction with a thermal bath with a mean 
photon number n, the atomic spontaneous emission with 
a decay rate 7 a for the upper level, and a loss of photons 
inside the fibers at a rate Kf. All dissipative effects can 
be described under the Markovian limit by standard Li- 
ouville superoperators so that the time evolution of the 
whole system density operator p(t) can be described by 
the following ME written in the Lindblad form: 



9 = 'I 



ri e ,P 



J=A,B,C 



(a) -r.(s)t 

c,J 



(2) 



where the non-Hermitian effective Hamiltonian is 



J 



(9)tA(3) 
c,J u c,J 



(3) 



The jump operators for the atoms are C 0; j = ^/tvtj, 
for the fiber losses Cf.j = ^/^Jfj, and for the cav- 
- \J K c (n + Tjcj (loss of a photon) and 



C^j — y/K c n,Cj (gain of a photon). The above ME can 



ity modes C). j 

be solved numerically by the Monte Carlo Wave Func- 
tion method [2^ |. From now on we consider dimension- 
less parameters, all scaled to the coupling constant gA, 
and times r = gAt- 

As a significative example, the implementation of our 
scheme may be realized in the optical regime of CQED 
by choosing a continuous variable (CV) entangled field 
for the subsystem (f) and two- level atoms as qubits (a). 
Each qubit is trapped in a one-sided optical cavity (c), 
where the radiation modes can be coupled to the cavity 
modes via optical fibers as in [30| . In optical cavities 
thermal noise is negligible (n = 0) , spontaneous emission 
can be effectively suppressed, and sin gle atoms can re- 
main trapped even for several seconds [24j . 

Here we focus on external field prepared in a qubit-like 
entangled state /5/(0) because this is the condition for 
high entanglement transfer for CV field [2l|. The gener- 
ation of photon number multimode entangled radiation 
was recently demonstrated [3lj . Under qubit-like behav- 
ior we can describe the entanglement of all three-qubit 
subsystems (a, c, f) by combining the information from 
tripartite negativity [32j |. entanglement witnesses [331 ] for 



the two inequivalent classes GHZ and W [3J|, and re- 
cently proposed criteria for separability (35|. In fact, the 
tripartite negativity (r) (a=a, c, f), defined as the 
geometric mean of the three bipartite negativities [36| . 
is an entanglement measure providing only a sufficient 
condition for entanglement detection, though its positiv- 
ity guarantees the GHZ-distillability that is an important 
feature in QI. 



III. STATE MAPPING AND TRIPARTITE 
ENTANGLEMENT TRANSFER FOR SINGLE 
MODE COUPLING 

A. Hamiltonian regime for external field in 
qubit-like pure states 

We first illustrate the Hamiltonian dynamics 
{{Hf, K c , 7a} <C 1) for the external field prepared 
in a qubit-like entangled pure state |^(0))/, atoms 
prepared in the lower state \ggg) a , and cavities in the 
vacuum state 1 000) c . By choosing vjk(t) — if J ^ K 
and vjj = gA we describe single-mode fibers each one 
propagating a mode of the entangled radiation. Wc 
show that under optimal conditions for mode matching 
it is possible to map |$(0))/ onto atomic and cavity 
mode states for suitable interaction times. Overall we 
are dealing with an interacting 9-qubit system, though 
the input field will be switched off at a time r // such 
that the atomic probability of excited state p e (r) reaches 
the maximum. Injected field switch-off can be obtained, 
e.g., by rotating fiber polarization. 

In Fig. [TJ we show numerical results for 
the external field prepared in the GHZ state 
|tf(0))/ = (|000)/ + |111> y)/ V / 2. In the time inter- 
val < t < r ff (transient regime) each flying qubit 
transfers its excitation to the cavity which in turn 
passes it onto the atom (see Fig. [TJ,). The cavity mode, 
simultaneously coupled to the external field and to the 
atom, exchanges energy according to a Tavis-Cummings 
dynamics at an effective frequency g^/2 [13, [33|- During 
the transient up to time r // = tt/V2 the mean photon 
number N^(t) = (c^c)(r) in each cavity completes a 
cycle. In the same period the atomic excitation proba- 
bility Peij) reaches its maximum value, while the input 
field has completely fed the cavity, i.e., its mean photon 
number N^\t) = (/V)( T ) vanishes. In Fig.[Tj3 we show 
that in the transient the atomic tripartite negativity is 
always positive and E^ a '(T ff) = 1, that is the value 
of the injected GHZ state. Until r // the dynamics 
maps the whole initial state |\P(0))/ ® |000) c ® \ggg)a 
onto the pure state |000)/ <g> |000) c <g> |#(0)) o , where 
|\&(0)) a is obtained from 1^(0))/ by the correspondence 
|0) / — > \g) a and |1)/ — > \e) a . This is confirmed in Fig. [TJ; 
by the time evolution of the purity /^"'(t) = Tr a [/5„(r)] 
and the fidelity F^(t) = a (*(0)|p a (r)|*(0)) a , where 
p a (r) is the atomic reduced density operator. As for 
the cavity mode dynamics we note that (see Fig. [T]3,d) 
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FIG. 1: Dynamics for the external field in a GHZ state: (a) 
iV (c) (dashed), iV (/) (dotted) and p e (solid); (b) E (a) for atoms 
(solid), cavity modes (dashed) and external field (dotted); (c) 
(solid) and F^ a) with <f> = (dashed), cf> = 7T (dotted); 



i (a) (solid) and F^ a) 
(d) ^ (c) (solid) and F 



(dotted). 



with 



-it/2 (dashed), 



+ 7T/2 



the local maximum of E^°'(T ff /2) does not correspond 
to a pure state, i.e. the initial state |^(0))/ cannot 
be mapped onto the cavity modes during the transient 
regime. The entanglement is only partially transferred 
to the cavity modes nevertheless allowing the building 
up of full atomic entanglement later on. This dynamics 
is quite different than in [2l| where the entangled field 
was mapped onto the cavity modes before the interaction 
with the atoms. 

At the end of the transient regime the external 
radiation is turned off and the subsequent dynamics 
is described by a triple JC ruled by oscillations at the 
vacuum Rabi frequency 2g, hence with a dimensionless 
period tt as shown by cavity mean photon number and 
atomic probability in Fig. QJ,. The purities /i( a,c )(r) 
in Figs. QJ;,d oscillate at a double frequency between 
pure entangled (maximum negativity) and separa- 
ble (zero negativity) states. In particular, at times 
T~m = Toff + m7r ( m = 0,1,2...) the atoms are in the 
entangled states (if |*(0)) o , where Of = (g) J e -^ & J &J 
is a local phase operator where 4> — Q {<p — tt) applies 
for even (odd) values of m, that are the peaks of 
E^(t) in Fig. [TJd. At times r„ = r o// + (n + ±)tt 
(n = 0,1,2...) the cavity mode states are obtained by 
applying = (g) J e- i * 6 J 6j , where <p = -f (+§) for 
even (odd) values of n, to the state |4 r (0)) c derived 
from 1^(0))/ by the correspondence |0)/ <-> |0) c and 
\l)f<-> |l) c . By choosing to turn off the external field at 
times shorter than r ff we find a progressive reduction 
in the entanglement transfer to the atomic and cavity 
subsystems, simulating the effect of non perfect cavity 
mirror transmittance. Up to 10% changes in the value of 
T ff, the fidelity F^(r ff) remains above 99.9%. We re- 



mark that the state mapping process can be obtained for 
any l^fO})/ written in a generalized Schmidt decompo- 
sition [3J] , as well as for mixed states as described below. 



B. Tripartite entanglement sudden death for 
external field in a Werner state 



For the injected field we consider the Werner state 
^(0) = (l-p)\GHZ){GHZ\ + p, {0<p< 1), because 
it is relevant in QI and it is possible to fully classify its 
entanglement as a function of parameter p. In fact, for 
< p < j the state belongs to the GHZ class and to 
the W class up to p = |. The tripartite negativity is 
positive up to 4/5, and due to the structure of /5/(0), 
the state is clearly inseparable under all bipartitions 
(INS), i.e. it cannot be written as convex combination 
of biseparable states. For 4/5 < p < 1 it is known that 
the state is fully separable 35|, [38| . 

The system dynamics can be divided into a transient 
and an oscillatory regime, and the state mapping of 
Pf(0) onto atoms (cavity modes) still occurs at times 
T m ( T n)- Out of these times the density matrices of all 
subsystems lose the form of a GHZ state mixed with 
white noise but still preserve invariance under all per- 
mutations of the three qubits and present only one non 
vanishing coherence as in /5/(0). This greatly helps us in 
the entanglement classifications in the plane (t, p) shown 
in Fig. [5] In fact, in the regions where E^ a \r) > but 
out of W class we can exclude the biseparability. The 
fully separability criteria in [3~5| are violated only where 
E^(t) > so that if E^(t) = the state may be fully 
separable or biseparable. Nevertheless, in the latter case 
the state should be symmetric and biseparable under 
all bipartitions and hence it is fully separable [39(. For 
any fixed value of p in the range < p < 4/5 we thus 
show the occurrence of entanglement sudden death and 
birth at the boundaries between fully separable and INS 
states. In particular, for < p < 4/7 we find genuine 
tripartite ESD and ESB phenomena. Note that, for 
a fixed value of p, the passage of atomic state during 
time evolution from W-class to GHZ-class (or viceversa) 
entangled states is permitted by the non-unitarity of 
the partial trace over non-atomic degrees of freedom (so 
that the overall operation on the initial three qubits is 
not SLOCC). We also notice that for times r > r e // we 
can solve exactly the triple JC dynamics, confirming our 
numerical results and generalizing (40| to mixed states. 

In Fig. [5Jd we see that, for increasing values of p, there 
is an increase of both the slope of B"(r) and the time 
interval of fully separability. In Fig. [2t we show in detail 
the transient dynamics of the tripartite negativities 
E^(t,p) (a = a,c,f) in the crucial region around 
T ff/2. We consider some values of p where the atoms 
exhibit in times different classes of entanglement. We 
see that for p = 0.2, where the input state has GHZ class 
entanglement, the ESB of subsystems (c),(a) anticipates 
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FIG. 2: ESD/ESB for external field in a GHZ state mixed 
with white noise, a) Regions in the plain (r,p) for atomic en- 
tanglement of type GHZ, W, INS, and fully separable (black), 
b) Sections E^ a \r) for selected values of p. c) Zoom on 
E (a \r f f /2) for field (dotted), cavity modes (dashed), and 
atoms (solid) with p — (black), p = 0.2 (blue), p — 0.4 
(green), p — 0.6 (red), d) Classification for cavity mode en- 
tanglement. 



the ESD of (f),(c), and there is an interval around r G ///2 
where all three subsystems are entangled (of INS-type). 
As p grows, hence the initial state becomes more noisy, 
the effects of ESD occur earlier and those of ESB later. 
For p = 0.4, involving W-class entanglement, only at 
most two subsystems are simultaneously entangled (first 
(f),(c) and then (c),(a)). For p = 0.6, involving only 
entanglement of INS-type, the cavity modes do not 
entangle at all (see Fig. [2Ji) ■ They physically mediate 
the discontinuous entanglement transfer from (f) to 
(a), where for p — > 4/5 the time interval without any 
entanglement increases while the entanglement level 
vanishes. 



C. Effect of dissipation on state mapping 

In the perspective of experimental implementation for 
QI purposes an important issue is the effect of dissipa- 
tion on both state mapping and entanglement transfer. 
For external field prepared in a GHZ pure state we first 
evaluated the effect of cavity decay rates in the range 
< k c < 0.5 for negligible values of all other decay rates. 
We consider as function of k c the behavior of the fideli- 
ties F^ a '(T min ) and the tripartite negativities E^ a ' (r mjn ) 
(a = c, /) at the first peaks (m = n = 0). In Fig. [3] we 
see that the above functions of k c can be well fitted by 
exponential functions, whose decay rates for the atomic 
subsystem are = 0.75, p£ — 1.09, and for the cavity 
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FIG. 3: Effect of cavity mode dissipation. At the first peaks 
T m ,n with m = n — we consider the tripartite negativities 
as function of k c = ft: F^(r ) (1), £ (a) (^o) (2), F< c »(t q ) 

(3),^(r„) (4). 

modes (3^) = 1.80, (3^ = 2.94. As expected, quantum 
state mapping and entanglement transfer are by far more 
efficient onto atomic than cavity qubits. For instance, if 
k c = 0.1 we obtain a state mapping onto the atoms (cav- 
ity modes) with a fidelity of S 0.93 (= 0.83). 

We can now add the further dissipative effect of atomic 
decay. For instance we find that, for an atomic decay rate 
7 a = 0.03 and in the presence of cavity decay with a rate 
k c = 0.1, the fidelity of the atomic (cavity mode) subsys- 
tem reduces by 4.4% (8.9%). 

Finally, we evaluate the effect of losses in the fibers 
used to inject the external field into each cavity. Clearly, 
this effect is relevant only up to the time r Q // = 2.22. 
We evaluated the effects of fiber decay rates kf up to 
1.0 for negligible values of atomic and cavity decay rates 
(k c « 1, 7 a << 1) (see Fig. 0]). We show the effect 
of parameter kf on cavity field mean photon number 
N^ c '(r ff /2) and atomic excitation probability p e (j ff) 
and we see that the amount of energy transferred to the 
atoms and to the cavity modes decreases exponentially 
for increasing values of kf, the decay rates are = 0.42 
and = 0.82, respectively. Also the behavior of the tripar- 
tite negativity E ( ^(t ) and fidelity F( q )(t ) at the first 
peak can be described versus kf by exponential functions 
whose decay rates are = 1.51 and = 1.95, respectively. 



IV. STATE MAPPING AND TRIPARTITE 
ENTANGLEMENT TRANSFER FOR 
MULTI-MODE COUPLING 

Finally, we consider multi-mode coupling of the 
external field to each cavity mode. For simplicity we 
choose equal coupling constants vj t K = 1/ j,k/9a ^ if 
K J and we consider values in the range 0—1.4. In 
the transient regime the dynamics is sharply modified 
with respect to the case of single mode fiber shown in 
Fig. [TJ By increasing the values of vj t K the period of 
energy exchange decreases from 27r/v / 2 to = 2.6. The 
maximum of cavity mode mean photon number grows 
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FIG. 4: Effect of fiber mode decay rate ft/ = Hf/gA for 
k c « 1, 7a << 1. We evaluate at time r Q // the atomic 
tripartite negativity E a (solid), the fidelity F a (dash-dot), and 
the atomic probability p e (dot), and at time r ///2 the cavity 
mean photon number (dash). 



up to N^ c ' = 0.41 whereas the maximum of atomic 
excitation probability decreases to p e = 0.24. The 
external field mean photon number does not vanish 
but it reaches a minimum, that can be always found 
between the two maxima of N^ c '(t) and p e (r), such 
that 0.002 < AT(/) < 0.02 changing i>j tK from 0.1 to 
1.4. We investigate the differences in the entanglement 
transfer for three selections of switching-off time r G // 
corresponding to the maximum of p e (r), the minimum 
of NW(t), and the maximum of N^(r). In Fig. Eh, 
we show the dependence of r /y on vj,k=£.j- Switching 
off the external field at times T Q ff corresponding to the 
maxima of p e (r), as in the previous case with single- 
mode fibers, we find (Fig. [SJd,c) that the maxima of 
tripartite negativities (r) after the transient regime 
reduce for increasing values of vj^k for both atomic and 
cavity mode subsystems. If we consider r a ff correspond- 
ing to the minimum of N^\t) (Fig. [SJi,e) we observe a 
small reduction of the peak values of £W(t), Finally, 
if we turn off the external field at the first maximum of 
the cavity field mean photon number we note that by 
increasing the values of vj y K it is possible to improve 
the entanglement transfer (Fig.[5fjg)- The peak value of 
tripartite negativity grows up to = 0.93 for vj^k = 1.4 
and the fidelity up to = 0.95 for both subsystems (a) and 
(c) . We remark that these values cannot be significantly 
increased for larger values of vj^k- In conclusion, for all 
the above choices of switching-off time r a ff we observe 
that, by increasing the values of vj t Ki the amount of 
entanglement that can be transferred to the cavity 
modes in the transient regime also increases. This is 
due to the fact that the amount of energy transferred to 
each cavity mode increases: in fact, the peak value of 
-/V( c )(t) progressively grows up from = 0.25 to = 0.41. 
Nevertheless, multimode coupling for larger values of 
T a ff results in a less favorite condition for entanglement 
transfer. 
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FIG. 5: Effect of multimode coupling, a) Dependence of 
T a ff on the coupling constants Vj t K for different choices of 
switching-off the external field: maximum of Pe{r) (o), min- 
imum of iV (/) (r) (x), and maximum of JV (c) (r) (+)• Tripar- 
tite negativities E^ a ' (a — a, c) for vjjx. = (solid gray), 
0.3 (dashed), 0.6 (dotted), 1.0 (dashed-dotted), and 1.4 (solid 
black): b,c) r // in the maximum of Pe(t); d,e) T ff in the 
minimum of JV'^(r); f,g) t // in the maximum of N( c \t). 



V. CONCLUSIONS 

In this paper we have addressed the transfer of quan- 
tum information and entanglement from a tripartite 
bosonic system to three localized qubits through their 
environments, also in the presence of external environ- 
ments. We considered an implementation in the optical 
regime of CQED based on CV photon number entan- 
gled fields and atomic qubits trapped in one-sided optical 
cavities, where the radiation modes can couple to cavity 
modes by optical fibers. 

In the nine-qubit transient regime the quantum state is 
mapped from tripartite entangled radiation to tripartite 
atomic system via the cavity modes. After the transient 
we switch off the external field. The subsequent triple JC 
dynamics, that we solved analytically (numerically) for 
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pure (mixed) input states, shows how the effect of map- 
ping can further affect the atom-cavity six qubits system. 
Its relevance is in the possible manipulation for QI pur- 
poses of entanglement stored in separate qubits of atomic 
or bosonic nature. Hence the interest to put quantitative 
limits dictated by cavity, atomic and fiber mode decays, 
that we evaluated at the times where the transfer proto- 
col is optimal. 

In the case of a GHZ input state mixed with white 
noise, we provide a full characterization of the separabil- 
ity properties of the tripartite subsystems. We can then 
show the occurrence of entanglement sudden death effects 
at the tripartite level, deriving the conditions for the re- 
peated occurrence of discontinuous exchange of quantum 
correlations among the tripartite subsystems. This is an 



issue of fundamental interest as well as worth investigat- 
ing for all applications in quantum information process- 
ing, remarkably computing and error correction, where 
disentanglement, which may be faster than decoherence, 
has to be carefully controlled. 

An extension and comparison to other types of entan- 
gled qubit-likc input fields and experimentally available 
CV fields will be presented elsewhere (42j . 
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